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ABSTRACT 

We examine the prospects for producing Nickel-56 from black hole accretion disks, by examin- 
ing a range of steady state disk models. We focus on relatively slowly accreting disks in the range 
of M = 0.05 Mq/s to M = 1 Mq/s, as are thought to be appropriate for the central engines 
of long-duration gamma-ray bursts. We find that significant amounts of Nickel-56 are produced 
over a wide range of parameter space. We discuss the influence of entropy, outflow timescale and 
initial disk position on mass fraction of Nickel-56 which is produced. We keep careful track of 
the weak interactions to ensure reliable calculations of the electron fraction, and discuss the role 
of the neutrinos. 

Subject headings: gamma ray:bursts-nucleosynthesis-accretion disks 



1. Introduction 

Over the last decade and a half there has been 
a dramatic increase in data available for studying 
long duration gamma-ray bursts (GRBs) . Besides 
the gamma rays, afterglows in the x-ray, opti- 
cal and radio are seen, and host galaxies have 
been identified. In a handful of cases "supernova 
bumps" are s een in the optical afterglows, e.g. 
GRB 980326 (jBloom et al.l ll999l) . Furthermore, 



spectra for Typ e Ic supernova have been identified 
in a few cases dGalama et al.l Il998t Irliorth et al ' 



2003; Kawab ata et all 2003: Mathcs on et al.ll200 . 
Stanek et al.l2003HMalesani et al.ll2004lThomsen et al 



2004) . These developments have lead to the theo- 
retical conjecture that all long duration gamma 
ray bursts are the result of the collapse of a 
massiv e star; for a review see Wooslev fc Blooml 
<|2006h . 



The temporal evolution of light curve spectra of 
traditional core collapse supernovae are driven by 
the decay of Nickel-56 to Cobalt-56 and the subse- 
quent decay of Cobalt-56 to Iron-56. The half lives 
for these decays are approximately six days and 
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seventy-seven days respectively. After these nu- 
clei beta decay, their daughter nuclei emit gamma 
rays which then thermalize in the ejecta. The en- 
ergy from these decays is thus re-emitted in the 
optical and observed as the supernova light curve. 
Thus the intensity of a supernova bump should 
be directly influenced by the amount of 56 Ni pro- 
duced in the object. On the order of tenths of 
solar masses of Nickel are needed. For example, 
it has been estimated that around 0.5M Q of 56 Ni 
is necessary to ex plain SN2003dh which acc ompa- 
nied GRB030329 dWooslev and Hegerll2003h . 

The commonly discussed model for the type 
of superno va that produces a GRB is the collap- 
sar m odel ( Wooslev 19931: MacFadven fe; Wooslev 
1999). In this case the star has too much rotation 



to produce a standard core collapse supernova. 
A hot ~ 10 MeV disk surrounding a black hole 
is produced at the center. Nickel-56 can be pro- 
duced in two different ways. It can be produced 
in explosive burning, when a shock wave propa- 
gates out through the sta r and pre-existing nuclei 
are then burned to 56 Ni (iMaeda fc Nomotoll2003t 



Fryer et al.|[200a iMaeda fc Tominagall2009h . This 



is the same mechanism as in traditional core col- 
lapse supernovae, but the elemental yields are dif- 
ferent since the explosion energy is larger. How- 
ever, Nickel-56 can also be produced in a disk 
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wind. This is a primary process where free neu- 
trons and protons are eje cted from the disk and 
combine to form nuclei rtSurman fe McLaughlin! 
20051: IPruet. Thompson, fc Hoffmanll2004l;ISurman et al 



hydrodynamics. In the inner regions of all of the 
disks, nuclei are dissociated into their constituent 
nucleons, and the ratio of neutrons to protons is 
subsequently set by the weak interactions 



20061) . 

Any assessment of 56 Ni production in GRBs 
must consider both mechanisms. The latter mech- 
anism is more difficult to model, since it re- 
quires a careful accounting of the neutrinos in 
the disk. In thi s paper we extend the prelim - 
inary studies o f Surman fc McLaughhnl (120051): 
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Pruet. Thompson, fc Hoffman ( 2 0041 ) : lPruet. Surman. fc 
( 2005 ): Surman et al.l (|2006l ). and examine the 



We describe the composition of the disks by the 
electron fraction Y e , where Y e = 1/(1 + n/p) and 
n and p are the number densities of neutrons and 
—pxataoSr-ifispectively. For disks with the lowest 
MsfeftuSfe^ates of those considered here, weak re- 



production of Nickel-56 from the outflow from 
a variety of disks. We discuss not only the im- 
portance of the electron fraction in determining 
the abundance yields, but also the role of en- 
tropy and outflow timesca le. Recent studies, e.g. 
Maeda fc Tominagal ( 2009! ) and references therein, 
have used estimates of ejected 56 Ni mass to con- 
strain the GRB central engine. This study pro- 
vides a reference for anyone considering Nickel-56 
production in models of long duration gamma-ray 
bursts. 

2. Disk and Outflow models 

In order to study the nucleosynthesis in the 
hot outflows from GRB accretion disks, we first 
need to determine the composition of and the 
neutrino fluxes emitted from the accretion disk 
itself. We therefore begin with existing hydro- 
dynamical black hole accretion disk (AD-BH) 
models with parameters appropriate for long- 
duration GRBs. For this study we have chosen the 
one-dimensional, steady-sta te AD-BH models of 
Chen fc Beloborodovl (|2007l) (CB). The CB mod- 
els are fully relativistic and incorporate important 
microphysics effects, such as neutrino trapping 
and an evolving neutron-to-proton ratio, neglected 
in other relativistic disk calculations. Here we in- 
vestigate a range of steady-state disk models with 
mass accretion rates 0.05 < m < 1, where fa = 1 
M Q /s. 

We calculate the initial nuclear composition of 
each d isk as described in Surman fc McLaughlin! 
([20041 ) . We find the disk compositions we calculate 
as a post-processing step m atch well with the disk 
compo sitions reported in IChen fc Beloborodov 
(|2007l ). which are calculated along with the disk 



action rates are slow even in the inner regions 
of the disks. In these cases the disk composition 
remains roughly balanced between neutrons and 
protons and Y e remains close to 0.5. For disks 
with higher accretion rates, conditions in the in- 
ner region of the disks are hotter and denser, and 
weak reactions are correspondingly faster. For 
disks with m < 1, electron capture on protons, 
the forward reaction of Eqn. 1, predominates, 
and the inner regions of the disks become quite 
neutron rich. Fig. [T] shows the electron fractions 
in the disk for two radial distances, ro = 50 km 
and 100 km, from the black hole for the disk 
models used in this work. A full discussion of 
the disk com positions of AD-BH can be found in , 
for example, Pruet. Wooslev. fc Hoffman! (2003 ): 



Fuiimoto et al.! (120031): ISurman fc McLaughlin! 



(|2004l ): IChen fc Beloborodovl (120071) . These disk 
electron fractions provide the starting points for 
the calculations of the nucleosynthesis from the 
disks. 

While calculating the composition of the disks, 
we additionally c alculate the neutrino fluxes from 
the disk as in ISurman fc McLaughlin! ([20041 ). 



Most of the AD-BHs considered in this work are 
optically thin to neutrinos and the emitted fluxes 
are primarily due to electron and positron cap- 
tures in the disk. Since the rate of electron cap- 
tures in the disks far exceeds the rate of positron 
captures for all disks with rh > 0.1, the electron 
neutrino flux emitted from the disk is significantly 
higher than the electron antineutrino flux. This 
has important consequences for the nucleosynthe- 
sis in the outflow. 

We follow the outflowin g mate rial as described 
in ISurman fc McLaughlin] (|2005l ). We take the 



outflow to be adiabatic and the velocity v of the 
outflowing material as a function of radial distance 
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Fig. 1. — Shows the electron fractions in the disk 
at radial distances of ro = 50 km (squares) and 
ro = 100 km (diamonds) for AD-BH with mass 
accretion rates rh. 

from the black hole r to be v = i>oo(l — r /r)P, 
with Vac = 10,000 km/s, r is the starting disk 
position, and f3 controls the rate of the acceler- 
ation of the material. We investigate a range of 
/3, 0.2 < f3 < 2.6, where smaller f3 corresponds 
to more rapid initial acceleration, and two disk 
starting positions, ro = 50 km and 100 km. We 
connect the evolution of the density to our veloc- 
ity prescription by setting the mass loss rate M 
of the outflow to be constant, first in cylindrical 
symmetry as the material lifts off vertically from 
the disk, and then in spherical symmetry as the 
material expands radially away at later times. We 
take the initial composition of the material to be 
that of the disk at the starting radius ro. We then 
follow the nuclear recombination u sing a full nu- 
clear network code as described in ISurman et al.l 
( 20061 ). A schematic of the outflow is shown in 
Fig.H 

While the entropy of the disk is quite low, sev- 
eral potential heating processes can increase the 
entropy of the material as it lifts off the disk. 
To make some estimate of what range in en- 
tropies is realistic for the wind, we consider the 
rate of change of entropy due to just one poten- 
tial heating source — neutrino captures in the out- 
flow. We start with the approximate expressions 
for the neutrino heatin g and cooling rates from 
Qian k. Wooslevl ( 19961 ) and calculate the change 
in entropy per baryon from energy conservation, 
recast for a neutrino-driven wind in Eqn (26) of 
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Fig. 2. — Shows a schematic of an outflow tra- 
jectory, labeled with the three stages of outflow 
nucleosynthesis. 



Qian fc Wooslevl (|l996l) as 
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where s is the entropy per baryon, tojv is the 
mass of a nucleon, T is the temperature, and q 
is the specific net heating rate due to neutrino 
and antineutrino interactions. When we apply 
this expression to our outflow velocity and den- 
sity profiles as described above, we find, for exam- 
ple, for an rh = 0.15 disk, the change in entropy 
per baryon in the outflow due to neutrino inter- 
actions ranges from 1 — 23fc, while for a rh = 1 
disk, the range is 4 — 58fc. We note that these 
estimates represent a lower limit, as additional 
processes (e.g., magnetic heating) are expected to 
play an important role. We therefore investigate 
outflow entropies up to somewhat larger values, 
10 < s/k < 100. 

3. Production of Nickel-56 

The element synthesis leading to the production 
of 56 Ni for a sample AD-BH outflow trajectory is 
depicted in Fig.[3l which shows the evolution of the 
mass fractions of neutrons, protons, alpha parti- 
cles, and 56 Ni as a function of decreasing temper- 
ature T 9 , where T 9 = T/10 9 K and T is the tem- 
perature in the outflow in Kelvins. In this exam- 
ple, the material is ejected from an AD-BH with 
rh = 0.5 at a starting radius of ro — 100 km; the 
initial composition is therefore quite neutron rich. 
As the material lifts off the disk, it transitions 
from the dense, low entropy, electron-degenerate 
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disk conditions to the higher entropy (s/k = 30 
in this example), hotter, accelerated outflow con- 
ditions. The higher entropy favors a balance of 
electrons and positrons such that positron cap- 
ture, the forward reaction of Eqn. 2, can compete 
with electron capture, and the electron fraction 
quickly rises. As the material moves away from 
the disk, it cools and expands, and weak reaction 
rates slow. Even so, the electron fraction contin- 
ues to evolve as the positron capture rate, and 
eventually the electron neutrino capture (the re- 
verse reaction of Eqn. 1) rate, are slightly faster 
than the electron capture rate. The positron cap- 
ture rate exceeds the electron capture rate because 
of the neutron-proton mass difference, and even- 
tually the electron neutrino capture rate exceeds 
both the electron and positron capture rates as the 
latter drop much more rapidly (~ T~ J ) above the 
disk. Therefore the material gradually becomes 
proton-rich prior to nuclear reassembly, as seen in 
Fig. |3 for T 9 > 10. This portion of the outflow 
trajectory, where the material is composed of free 
nucleons and their ratio is adjusted by weak inter- 
actions, is labeled (a) in Figs. [2] and El 

In the first stage of nuclear reassembly, labeled 
(b) in Figs. [5] and [31 protons and neutrons combine 
into alpha particles ( 4 He nuclei) . Conversion into 
alphas is most efficient at Y e = 0.5, where equal 
numbers of protons and neutrons are present. In 
this example, the material is slightly proton-rich, 
so all of the neutrons and a matching number of 
protons combine into alphas, and a small fraction 
of excess protons are left over. In the second stage 
of nuclear reassembly, labeled (c) in Figs. [2] and [3j 
alpha particles combine via the triple alpha pro- 
cess to produce 12 C, and a series of alpha captures 
on carbon culminate in the production of 56 Ni. 
The triple-alpha process and the subsequent alpha 
captures are highly temperature-dependent, and 
so the efficiency of the conversion of alphas to 56 Ni 
depends strongly on the evolution of the thermo- 
dynamic conditions in the outflow. If the outflow 
is coasting at this point such that the tempera- 
ture drops slowly, most of the alphas will convert 
to nickel; in the example considered here, however, 
this isn't quite the case, and a bit more than one 
half of the alphas end up in 56 Ni. 

The production of 56 Ni can therefore depend on 
a number of parameters: the mass accretion rate 
of the AD-BH, the point in the AD-BH at which 
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Fig. 3. — Shows the mass fractions of neutrons 
(solid line), protons (dashed line), alphas (dotted 
line) and 56 Ni (dot-dashed line) as a function of 
decreasing temperature Tg for an outflow trajec- 
tory with parameters r$ — 100 km, s/k = 30, 
and /3 = 0.8 from an AD-BH with ra = 0.5. The 
shaded areas refer to the portions of the outflow 
trajectory where (a) the ratio of free protons to 
neutrons is set by weak interactions (dark shaded 
region) , (b) protons and neutrons combine to form 
alphas (light shaded region), and (c) alphas com- 
bine to make 56 Ni (unshaded region). 

the outflow originates, the entropy per baryon and 
the acceleration of the outflow, and the neutrino 
interactions in the outflow. The influence of each 
of these parameters in turn is described in sections 
l3~TH3~3l below. 

3.1. Variation with entropy 

The production of nickel depends sensitively on 
the entropy per baryon in the outflowing material. 
Fig. 0] shows the final mass fraction of 56 Ni in the 
disk outflow as a function of the mass accretion 
rate m of the AD-BH and the entropy s/k in the 
outflow, for a sample set of trajectories with start- 
ing disk position r = 100 km and wind parameter 
f3 = 0.2. We see vigorous production of 56 Ni over 
the majority of the parameter space. Interestingly, 
for m > 0.2, the production of nickel does not ap- 
pear to depend strongly on the accretion rate of 
the AD-BH and is instead a much stronger func- 
tion of the entropy per baryon in the outflow. 
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Fig. 4. — Shows nickel production as a function of 
the entropy s/k of the outflowing material and the 
mass accretion rate rh of the disk from which the 
material originates. All trajectories shown have a 
starting disk position of tq = 100 km and wind 
parameter ft = 0.2, which corresponds to rapid 
initial acceleration. The shaded regions indicate 
mass fractions of 56 Ni in excess of 0.5 (darkest 
shaded region), 0.4, 0.3, 0.2, 0.1, 0.05, and less 
than 0.05 (unshaded region). 

3.1.1. m>0.2 

First we consider the nickel production in out- 
flows from the AD-BHs with rh > 0.2. In the out- 
flow trajectories with s/k < 20, the entropy is not 
much higher than in the disk itself and so the ma- 
terial retains some of the neutron richness of the 
disk. In these cases, the neutron excess remain- 
ing after alpha-particle assembly leads to neutron 
captures on the alpha-capture chain nuclei, and 
heavier isotopes of nickel predominate. Outflow 
trajectories with s/k > 20 are sufficiently high in 
entropy such that the neutron richness of the disk 
is erased by positron capture in the outflow. The 
most efficient production of 56 Ni occurs just at this 
transition between neutron-rich and proton-rich, 
where Y e = 0.5; for the trajectories considered in 
Fig. [4] this occurs when s/k ~ 30. At even higher 
entropies, the assembly into alphas is shifted to 
later times, and the weak reactions on free nucle- 
ons have more time to drive the material proton 
rich. Therefore as we consider outflow trajectories 
with various entropies, we find the higher entropy 
cases produce higher Y e and thus less efficient pro- 
duction of 56 Ni. 



3.1.2. m < 0.2 

The outflows from AD-BHs with rh < 0.2 dif- 
fer slightly from the above picture. Disks with 
lower mass accretion rates have lower central disk 
temperatures and densities. The most important 
effect of this for disks with 0.1 < rh < 0.2 is that 
the outflows have lower starting temperatures and 
densities. This means the weak rates are not nec- 
essarily fast enough for the electron fraction to 
reequilibrate in the outflow prior to nuclear re- 
assembly, and so the neutron richness of the disk 
is retained to higher entropies. For the lowest 
mass accretion rate disks considered in this work, 
0.05 < rh < 0.1, the central disk temperatures and 
densities are sufficiently low that the composition 
of the disk remains fairly balanced between neu- 
trons and protons, as shown in Fig. [TJ The out- 
flows from these disks therefore have electron frac- 
tions close to a half — and correspondingly strong 
56 Ni production — at lower entropies. 

3.2. Variation with ft 

While the entropy per baryon in the outflow 
plays a central role in setting the electron fraction, 
another important ingredient is the time available 
for each stage of the nucleosynthesis. In our out- 
flow parameterization this is set by the wind pa- 
rameter ft, which determines the initial accelera- 
tion of the outflow material. Fig. [5J similarly to 
Fig. HI shows the final mass fraction of nickel in 
the outflowing material as a function of the mass 
accretion rate rh of the AD-BH and the entropy 
s/k in the outflow, for a sample set of trajectories 
with starting disk position tq — 100 km and wind 
parameters ft = 0.8, 1.4, 2.0, and 2.6. As shown in 
Fig.0 increases in ft decrease the overall efficiency 
of 56 Ni production and, for high values of ft, shift 
the entropy per baryon s/k of the most efficient 
production. 

In order to understand the decrease in efficiency 
with increasing ft, consider that the most robust 
56 Ni production occurs when the outflow condi- 
tions evolve (1) sufficiently rapidly at early times 
to prevent the material from becoming too proton- 
rich and (2) sufficiently slowly at late times to fa- 
cilitate the complete conversion of alphas to 56 Ni. 
Small values of ft correspond to rapid initial accel- 
eration where the material quickly attains coast- 
ing speed, while outflows with large values of ft 
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Fig. 5. — Similar to Fig. [4j but for wind parame- 
ters = 0.8, 1.4, 2.0, and 2.6. 

accelerate slowly at first and evolve more rapidly 
at later times. This effect is illustrated in Table 
1, which shows the average values of the outflow 
velocity r divided by the radial distance of the ma- 
terial r for two key points in the evolution of the 
outflow — at early times when the electron fraction 
is set by weak interactions and late times during 
nuclear reassembly. The value of f/r serves as an 
indicator of how rapidly conditions evolve during 
these key periods; robust 56 Ni production there- 
fore requires large values of f/r at early times and 
small values at late times. As Table 1 shows, these 
conditions are satisified in trajectories with the 
smallest values of j3. Trajectories with larger /3 
evolve too rapidly at late times to efficiently syn- 
thesize nickel, and, as shown in the table, most of 
the material is left as alpha particles. 

For very high j3, f3 > 2.0, the most efficient 56 Ni 
production shifts to lower entropies. This effect is 
due to the influence of neutrino interactions in the 
outflow. As shown in Table 1, the evolution of 
the outflow for high (3 is slow at early times, when 



the material is close to the disk. The high /3 tra- 
jectories, therefore, receive the highest neutrino 
flucncc. Since the electron neutrino flux from the 
disk is significantly higher than the electron an- 
tincutrino flux, the net influence of the neutrinos 
is to drive the material proton rich via the reverse 
reaction in Eqn. 1. The low entropy trajectories 
that, in the absence of the neutrino interactions, 
would be neutron rich become proton rich and can 
therefore effectively produce 56 Ni, as illustrated in 
Fig.l 
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Fig. 6. — The top panel is identical to the third 
panel of Fig. [5l which shows nickel production as 
in Fig.|4]for wind parameters /? = 2.0 and ro = 100 
km. The bottom panel is identical except the neu- 
trino interactions in the outflow have been turned 
off. 



3.3. Variation with starting disk position 

While the outflow nucleosynthesis is primarily 
dependent on the outflow conditions as described 
above, the initial conditions — the local disk con- 
ditions at the point of origin of the outflow — also 
have some bearing on the final nickel production. 
Fig. similarly to Fig. 21 shows the mass fraction 
of nickel in the outflowing material as a function 
of the mass accretion rate m of the AD-BH and 
the entropy s/k in the outflow, for a sample set 
of trajectories with /3 — 1.4 and two starting disk 
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positions, tq = 50 km and 100 km. Nickel produc- 
tion is similar in both cases, though trajectories 
originating from the inner disk region are some- 
what less efficient. This is due to two separate 
effects — the influence of neutrino interactions and 
the starting temperature of the material. In the 
disk itself, the temperatures and densities increase 
with decreasing radius, and correspondingly the 
rates of weak reactions, and therefore the emis- 
sion of neutrinos, also increase. Thus the material 
outflowing from smaller radii will have higher ini- 
tial temperatures and densities than outflows from 
larger disk radii, and the neutrino flux the mate- 
rial is exposed to in the outflow will be higher. The 
latter effect tends to increase the electron fraction 
of the material; for trajectories that are already 
proton rich this decreases the efficiency of alpha 
particle formation. The former effect — the larger 
of the two — postpones alpha particle reassembly 
to later times and farther distances in the out- 
flow. As shown in Table 1, the velocity divided by 
radial distance r/r is faster during the epoch of al- 
pha particle reassembly for smaller starting radii, 
resulting in less efficient nickel production. 
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Fig. 7. — Similar to Fig- SJ but for wind parameter 
(3 = 1.4 and starting disk positions ro = 100 km 
and 50 km. 



4. Additional nucleosynthesis 

While in this work our focus is on the produc- 
tion of Nickel-56, we note that other interesting 
nuclear species may also be synthesized in the out- 
flow. Trajectories with strong Nickel-56 produc- 
tion can make a number of intermediate mass nu- 
clei that are underproduced in galactic chemical 
evolution studies, such as 45 Sc, 49 T i, and 64 Zn 
(jPruet. Surman. fc McLaughlin! 120051 ) . Low en- 



tropy (s/k < 20) trajectories from disks with m > 
0.1 are moderately neutron rich; these outflows 
produce heavier nickel isotopes and N — 50 nu- 
clei, and can potentially make some light r-process 
nuclei. High entropy (s/k > 50) trajectories with 
slow initial acceleration can make some interesting 
light p-process nuclei such as 74 Se, 78 Kr, 84 Sr, and 
92 Mo. 

An example of the latter two types of nucle- 
osynthesis is shown in Fig. [FJ which compares the 
final mass fraction of the light p-nucleus 92 Mo pro- 
duced in the outflow to its solar mass fraction. The 
origin of this nucleus has long been a mystery, 
thou gh studies of the newly-discovered vy pro- 



cess (jFrohlich et all 1200a iPruet et all 120061 ) sug 



gest early supernova neutrino-driven winds to be 
a strong candidate site for its production. Here we 
examine its production in the context of our GRB 
outflow models. 

In the case of the moderately accelerating out- 
flow, j3 = 1.4, we see strong 92 Mo production 
for low entropies, where the material is moder- 
ately neutron-rich. However, here the synthesis of 
92 Mo is accompanied by excessive overproduction 
of some N = 50 closed shell nuclei such as 88 Sr, 
89 Y, and 90 Zr. The overproduction of N = 50 nu- 
clei is so large (greater than a factor of 10 5 over 
solar in some cases) for these trajectories that if 
one knew precisely the amount of disk outflow and 
the rate at which accretion disks are formed, one 
could in principle constrain the disk wind condi- 
tions from this information. 

Alternately, we see that for the case of the out- 
flows with slow initial acceleration, (3 = 2.6, the 
greatest 92 Mo production occurs at high entropies, 
from disks with the most vigorous neutrino emis- 
sion. These trajectories have electron fractions 
greater than 0.6 and experience the highest neu- 
trino fluence, both requirements for a successful vp 
process. The production of 92 Mo and other light 
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Fig. 8. — Shows the final mass fraction of light 
p-process nucleus 92 Mo divided by its solar value, 
for wind parameters /3 — 1.4 (bottom panel) and 
(3 = 2.6 (top panel), and tq = 100 km. The shaded 
regions indicate 92 Mo production ratios in excess 
of 10000 (darkest shaded region), 1000, 100, and 
10 (lightest shaded region). 



p nuclei via a vp process in GRB disk outflows has 

(l2010h . 



been discussed inlKizivat et al 



5. Alternate outflow model 

Finally, for comparison we consider an alter- 
nate, spherical outflow model in which the evolu- 
tion of the electron fraction is tied more closely 
to the hydrodynamics. We construct the model 
by first constraining the outflow to follow ra- 
dial streamlines with zero vorticity. This, along 
with the steady flow approximation, implies the 
Bernoulli function is constant along the streamline 
that defines the outflow trajectory. In its simplest 
form, the Bernoulli function is 



B = -v 2 



h + $, 



(4) 



where v is the velocity, h is the specific enthalpy, 
and $ is the gravitational potential due to the 
black hole, —GMbh/t. We write the specific en- 



^ riifii + (kT)(s/k)n B 



(5) 



where is the baryon number density, and 
and ^ are the number density and chemical po- 
tential, respectively, of species i, where the sum 
over species includes electrons, positrons, neu- 
trons, protons, alpha particles, and an average 
heavy nucleus. We solve for the outflow hydro- 
dynamics by further assuming adiabatic flow and 
a constant mass outflow rate, M = 4nr 2 pv. We 
evolve the composition of the outflowing material 
as described in Sec. 2. 

The spherical outflow model differs from the 
parameterized wind primarily in the shape of the 
outflow trajectory — which is purely radial with no 
initial vertical piece — and the significantly slower 
acceleration of the outflowing material. Both of 
these effects impact the extent to which neutrino 
interactions influence the resulting nucleosynthe- 
sis, though in opposite directions. A purely radial 
trajectory moves material away from the central 
disk region, which is the source of the hottest neu- 
trinos, more directly than initially-vertical flow. 
However, this effect is dwarfed here by the effect 
of the dramatically larger neutrino fluence that 
arises due to the slower acceleration. The ini- 
tial flow is so much slower that the electron frac- 
tion re-equilibrates in the outflow in every case, 
to values above 0.5. For the low mass accretion 
rate disks, this results in vigorous 56 Ni produc- 
tion, as shown in Fig. [9] However, in higher accre- 
tion rate disks that have non-negligible antineu- 
trino emission, the production of nickel shifts to 
heavier isotopes. This is due to a mini-z/p-likc 
process, whereby antineutrino interactions on the 
free protons remaining after the onset of heavy el- 
ement formation create neutrons that are quickly 
captured by nuclei, primarily the dominant 56 Ni. 
The result is a shift from 56 Ni to 58 Ni and finally 
to 60 Ni with increasing disk m, as shown in Fig.[9j 
compared to the almost constant production of 
56 Ni that results when neutrino interactions are 
neglected. 

6. Conclusion 

Nickel production in the outflows from black 
hole accretion disks may be an important com- 
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Fig. 9. — Shows the final mass fractions of 56 Ni 
(solid red line), 58 Ni (long-dashed green line), and 
60 Ni (short-dashed blue line) as a function of disk 
mass accretion rate to, for the spherical outflow 
model with s/k = 40 and ro = 50 km, compared 
to the final mass fraction of 56 Ni (dotted red line) 
obtained with the same outflow model but without 
neutrino interactions. 

ponent of the total GRB Nickel-56 production. 
We have outlined the ways in which 56 Ni can be 
formed by black hole accretion disks under a va- 
riety of conditions, taking careful account of the 
neutrino interactions. 

We find that entropy is a strong diagnostic 
across a broad range of accretion disk models. 
Nickel production appears to be optimal for only 
modest increases in entropy over the in-disk val- 
ues; for the disks studied here with m > 0.1, this 
can be accounted for via neutrino heating alone. 
We note that nickel can also be produced via ther- 
monuclear burning from a shock that goes through 
a star, if the disk is associated with a hypernova 
explosion. However, in models without this addi- 
tional production mechanism, one could constrain 
the required conditions for nickel production in 
the object. For example, if accretion disk outflows 
are to be a significant fraction of the overall GRB 
nickel production, the contribution of other pro- 
cesses, such as magnetic heating, to the entropy 
increase in the outflow from to > 0.1 disks should 
not greatly exceed that of the neutrino interac- 
tions. 

However, conclusions about the mass fraction 
of Nickel produced cannot be drawn from the en- 
tropy alone. Nickel formation is optimized when 
the outflow timescale, in particular the ratio oir/r 
in the outflow, is fast at the beginning stages of nu- 



cleosynthesis and slow later. Furthermore, in the 
models considered here, outflows from an inter- 
mediate disk radius, ro ~ 100 km, more efficiently 
produced nickel than material ejected closer to the 
black hole. Under the optimal conditions up to 
around 50% of the material ejected from black hole 
accretion disks will become Nickel-56. These re- 
sults, when coupled with future wind model cal- 
culations and observations, can ultimately be used 
to constrain the outflow conditions of GRB disks. 
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r = 100 km 


p 


f/r, (a) 


f/r, (c) 


A"56Nj 


X a 


0.2 


105 


2.3 


0.612 


0.314 


0.8 


32 


6.1 


0.488 


0.449 


1.4 


4.8 


14 


0.327 


0.590 


2.0 


0.62 


28 


0.189 


0.667 


2.6 


0.07 


32 


0.131 


0.635 




ro = 50 km 


(3 


f/r, (a) 


f/r, (c) 




X a 


0.2 


120 


3.1 


0.582 


0.355 


0.8 


52 


9.6 


0.422 


0.523 


1.4 


5.2 


27 


0.229 


0.689 


2.0 


0.62 


58 


0.079 


0.752 


2.6 


0.06 


64 


0.048 


0.705 



Table 1: Average values of velocity divided by ra- 
dial distance f/r in units of s _1 for two sets of 
trajectories with entropy s/k = 30 and starting 
disk positions ro = 100 km and ro = 50 km from 
an AD-BH with m = 0.5. Each set of trajecto- 
ries includes five values of wind parameter /3, and 
the average f/r is calculated twice: (a) at early 
times as the electron fraction is set (Tg > 10) 
and (c) at later times during nuclear reassembly 
(6 > Tg > 2), where (a) and (c) correspond to the 
portions of the trajectory depicted in Fig. [2] Also 
included are the final mass fractions of 56 Ni and 
alpha particles for each trajectory. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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